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We analyze the relaxation of a superconducting flux qubit during measurement. The qubit state 
is measured with a nonlinear oscillator driven across the threshold of bifurcation, acting as a switch- 
ing dispersive detector. This readout scheme is of quantum non-demolition type. Two successive 
readouts are used to analyze the evolution of the qubit and the detector during the measurement. 
We introduce a simple transition rate model to characterize the qubit relaxation and the detec- 
tor switching process. Corrected for qubit relaxation the readout fidelity is at least 95%. Qubit 
relaxation strongly depends on the driving strength and the state of the oscillator. 

PACS numbers: 03.67.Lx, 85.25. Cp , 85.25.Dq 



Superconducting qubits are quantum systems based 
on microfabricated superconducting circuits with one or 
more Josephson junctions as nonlinear elements [T]. They 
are artificial quantum systems, with properties that can 
be defined by design of the mesoscopic parameters of 
the circuit. In superconducting qubits quantum state 
readout is of considerable interest, since the fabricated 
nature of qubit and detector allows full control of the 
qubit-detector coupling strength. Consequently, aspects 
of quantum measurement can be experimentally inves- 
tigated that are commonly not easy to access. These 
include the realization of high fidelity [2] and projective 
[3J 15 measurements, partial measurements [5], and the 
continuous observation of qubit dynamics [HI |7] . In addi- 
tion, state readout is a subject relevant to quantum com- 
puting: projective measurements are an essential part of 
protocols for quantum information processing. 

Qubit state readout can be performed in various ways. 
In dispersive readout the qubit is coupled to an oscillator, 
with a quadratic type of interaction. As a result of this 
nonlinear coupling, the resonance frequency of the oscil- 
lator becomes qubit-state dependent |7J[^. The state of 
the qubit can thus be inferred from a measurement of the 
properties of the oscillator. 

Nonlinear switching detectors are very attractive as 
they are able to amplify the information extracted from 
the qubit, leading to very fast readout with high fi- 
delity. Here we present a detailed experimental analysis 
of switching dispersive readout of a superconducting fiux 
qubit with a nonlinear oscillator. We introduce a simple 
model that allows to characterize the detector switching 
process and the qubit relaxation induced by the oper- 
ation of the detector. We find that the main source of 
measurement error is qubit relaxation induced by the op- 
eration of the detector. 

The persistent current flux qubit circuit [9] is a super- 
conducting loop interrupted by three Josephson junctions 
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FIG. 1: (a) Qubit and readout circuit, (b) Bistability dia- 
gram of the oscillator for the two qubit states |gi > or |e >. (c) 
Oscillator driving amplitude for qubit readout, (d) Oscillator 
switching probability, P{H), as a function of the switching 
plateau driving amplitude Pw for the qubit prepared initially 
with a Rabi pulse, 0{g), n{e) , 37r/2(s). (e) First derivative 
ofP{H){P^). 



(see figure |l|i) . Biased with an external flux close to 
half a flux quantum <&o/2, it behaves as a quantum two- 
level system. The ground state \g > and the excited 
state |e > are quantum superpositions of two oppositely 
circulating persistent currents (±/p). In the basis of the 
current operator / = Ip&z, the Hamiltonian of the flux 
qubit is: H = i(ea^ -I- Aa^), where e = 2Ip{^qb - l^o) 
and A is the quantum tunneling energy between the two 
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FIG. 2: (a) Oscillator driving sequence used to analyze the 
qubit relaxation. In the first driving pulse the oscillator is set 
either in L ox H (see text), (b) Decay of the probability for 
the measurement outcome L (due to the qubit relaxation), (c) 
Qubit relaxation rates versus the oscillator driving amplitude 
Idrive for the osclUator in the state L or H. 



current states. Here A/h = 5 GHz and the qubit is op- 
erated at a frequency Fqb — + A^/Zi = 14.2 GHz. 

Our detector is an oscillator formed by the inductance 
of a DC-SQUID and a capacitor. The inductance of the 
DC-SQUID and consequently the resonance frequency 
of the oscillator Fq depend on the flux enclosed in 
the SQUID loop. The oscillator is operated at a fre- 
quency Ff) — 1.5 GHz and has a quality factor Q = 20. 
The SQUID inductance is nonlinearly dependent on the 
SQUID current. Due to this nonlinear! ty, when the os- 
cillator is driven at a frequency F < Fq{1 ~ it can 
switch between a state of low oscillation amplitude (la- 
beled L) and a state of high oscillation amplitude (la- 
beled H) [TO]. Three situations are possible depending 
on the amplitude of the driving current lac compared to 
the two bifurcation currents Ibi{F) < Ii,2{F). For weak 
driving lac < hi the oscillator is in the L state. For 
strong driving lac > Ib2 the oscillator is in the H state. 
For intermediate driving, < lac < Ib2 the oscillator is 
bistable and can be in either of the L, H states. 

Due to the flux-dependence of the SQUID inductance, 
the upper bifurcation current Ib2 is highly sensitive to the 
flux (^sq- As the expectation value of the flux generated 
by the qubit is different for the two energy eigenstates, 
Ib2 depends on the qubit state. In particular we choose 
the qubit bias such that Ib2{g) < ^62 (e). For the readout 
of the flux qubit, the oscillator driving amplitude Isw is 
increased to a value Ib2{9) < Isw < h2{e) (see flgures[l]3 
and[l|;), such that the oscillator switches to the H state 
with a high probability if the qubit is in \g >, while it 
stays in the L state if the qubit is |e >. This first time 
interval (the switching plateau) with duration t^w consti- 
tutes the actual measurement interaction. To optimally 
discriminate between the two oscillator states L and H, 
noise from the detection electronics needs to be averaged 



out. This is performed during the holding plateau (thoid) 
with the driving amphtude decreased such that both os- 
cillator states can be maintained without switching or 
retrapping. 

It should be noted that the large difference between 
Fqb and Fq inhibits energy exchange between the qubit 
and the oscillator. The states |e > and \g > are thus 
preserved during the measurement. 

As shown in figures and [l^, the two bifurcation 
currents Ib2{g) and /b2(e) can be remarkably well re- 
solved. The qubit readout is performed at an ampli- 
tude 1°^ = 115 nA, where the switching probability is 
pa{H) = 99.7% if the qubit is in \g > and P^(i?) = 
14.6% if the qubit is in |e >, resulting in a readout con- 
trast of 85%. The main loss of readout contrast occurs 
when the qubit is in |e >, suggesting that the readout fi- 
delity is limited by the qubit relaxation before or during 
the measurement. 

In practical flux qubits one finds a very irregular de- 
pendence of the relaxation on flux bias due to electromag- 
netic modes and 'natural' two-level systems. The qubit 
is operated at a flux bias where the relaxation rate is lo- 
cally minimal. During measurement, the SQUID trans- 
port current varies in time. By second-order processes, 
this oscillation shifts the average value of the circulating 
current thus shifting the qubit bias. At that new point, 
relaxation is likely to be faster. More seriously, the qubit 
flux bias is swept at the oscillator driving frequency. In 
the H state, the flux is modulated over a range as large as 
5 m<i>o corresponding to a sweep of the energy splitting 
over 5 GHz. The qubit can thus be swept through re- 
gions where the relaxation rate is much higher. Moreover, 
when the oscillator is driven into its nonlinear regime new 
channels of relaxation might open where qubit energy is 
transferred directly to the oscillator [12]. 

We first measure the qubit relaxation under conditions 
where the oscillator is fixed in either the L or the H state. 
The qubit is initially prepared in |e >. Next the oscillator 
is set either in the H state or in the L state. To prepare 
H a short high driving pulse lac > Ib2 (dashed line in 
figure |2^) is applied. Subsequently the oscillator is driven 
for a time t drive at an amplitude Idrive- Afterwards the 
qubit state is read out with a regular measurement pulse. 
Figure as an example gives the decay in time of the 
probability for readout in the L state when the oscillator 
is driven with an amplitude Idrive = Ihoid- This decay 
is exponential and is due to qubit relaxation. The two 
qubit relaxation rates and F^^ for the two states of the 
oscillator are significantly different. Figure [2]; shows the 
dependence of the qubit relaxation rates on the oscillator 
driving amplitude for Idrive < hold and hrive > hold 
when the oscillator is in the L or H state, respectively. 
Whereas F^ is almost constant, F^^ increases with the 
driving amplitude. 

Similarly, we characterize the qubit relaxation during 
the rising part of the readout pulse with the effective re- 



3 



laxation rate T"'"' obtained from an exponential decay 
fitting of P{L) as a function of the rise time Uise- Cor- 
rected for the relaxation during the rising part of readout 
pulse (5%) and initial qubit preparation errors (5%) [T3], 
the readout fidelity is / = 95%. The remaining errors 
occur during the switching plateau. 

The approach used so far is well suited to analyze the 
qubit relaxation for a driving amplitude where the os- 
cillator is in a stable state. However for higher driving 
amplitude Idrive > Ihoid and especially at the driving 
amplitude during the switching plateau, the L state is 
metastable and can switch to the H state. Therefore 
to analyze the qubit relaxation, the oscillator switching 
process needs to be included. 

In the following we analyze the oscillator switching 
and the qubit relaxation during the switching plateau. 
The solid line in figure [SJd shows the oscillator switch- 
ing probability when the qubit is in jt; > as a function 
of the duration of the switching plateau tsw We distin- 
guish two regimes, indicated as I and II. The boundary of 
the two regimes is after a time of about 10 ns, indicated 
as tj/jj. In regime I the switching probability increases 
very fast. At the origin of the switching plateau, the 
oscillator driving amplitude is increased in about 1 ns, 
probably leading to non-adiabatic effects. In regime II 
the switching probability increases with a constant rate. 
The switching rate Tgw from the L state to the H state 
depends on the height of the effective potential barrier 
[m A[/ = Uo[l - (/ac//fc2)^]^/^ between L and H. For a 
driving current lac close to the upper bifurcation current 
If,2, the oscillator switching rate Tsw increases strongly 
by a few orders of magnitude. As Ib2 is different for 
the two qubit states, T^w strongly depends on the qubit 
state, which is the principle of the measurement. As 
42(5) < h2{e), it follows that Tf^ « TJ^ (figure [s]:) . 

We describe the qubit and the oscillator as a 4-state 
system {g,L), {g,H), {e,L) and {e,H), with 4 cor- 
responding occupation probabilities. Due to normal- 
ization, only 3 probabilities are independent. We use 
P{H) = P{e,H) + P{g,H), P{e,H) and P{e,L). In 
regime II, oscillator switching {L —> H) and qubit relax- 
ation {e ^ g) are described with the set of four rates 
rf^,r^^ and Tf,Tf (figure [sji) . Oscillator retrapping 
{H — !■ L) and qubit excitation (g — > e) are negligible. 
The evolution of the occupation probabilities is given by: 



dP{e,L) 
dt 

dP(e,H) 
dt 

dP{g,L) 
dt 

dP{g,H) 
dt 



-P{e,Lm^,+V\) 
Pie,L)Tt^-P{e,H)Tf 
Pie,L)Tf-P{g,L)TI^ 
P{g,L)TI^ + P{e,H)Tf 



(1) 



For a given initial qubit state i (g or e) at the ori- 
gin of the switching plateau t = 0, we denote Pq{H), 
PQ{e,H), and PQ{e,L) as the initial conditions for the 
set of equations (jlj . They quantify the evolution of the 
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FIG. 3: (a) Oscillator driving amplitude, (b) Switching 
probability as a function of the duration of the switching 
plateau (solid line), (c) Effective potential barrier between 
L and H, for the qubit state \g > or |e >. (d) Schematic 
of the model used to describe the evolution of the qubit and 
the oscillator during the readout pulse, including qubit re- 
laxation rates ,r]^ and oscillator switching rates rj^jFf^. 
Pg{H),PSiH),PS{e,H),PS{e,L) describe the regime I (see 
text). 



qubit and the oscillator evolution during regime I. We 
choose to define the initial conditions at i = 0, and not 
at t — tj/jj. They are obtained by extrapolating P{H), 
P{e, H) and P{g, L) in regime II to t = 0. In the more 
general case of an initial qubit state with an occupation 
probability of |e >: P{e) = x, the initial conditions 
are: Po{H) ^ xP§{H) + {I - x)P^{H) (and similarly 
for Po{e,H) and Pole.L)). 

If the qubit is initially in \g >, the only relevant rate is 
Ff^ and the set of equations ([ij reduces to: dP{H) /dt — 
[1 — P(iJ)]Ff^. For a given switching plateau driving 
amplitude Isw, Po{H) and Ff^ are extracted by fitting 
the switching probability P{H) as a function of with 
P{H) = !-[!- P^{H)] exp (-^siorf^) (as shown by the 
dashed line in figure |3|d) . 

The driving amplitude 1°^ used for the measurement 
is higher than Ib2{g)- In that case the L state does not 
exist when the qubit is in |g >, so qubit relaxation di- 
rectly causes oscillator switching. The state (g, L) is 
eliminated from the set of equations ([T]) by assuming 
that FJ" corresponds to the transition (e, L) {g, H) 
shown by the dashed arrow in figure [SJi. If the qubit is 
in |e >, two switching processes are possible, either due 
to qubit relaxation (rate F^"), or to switching while the 
qubit is excited (rate F^^). The sum F^^ = T^ + T1^ and 
Pq{H) can be extracted from P{H), given by: P{H) = 
1 — [1 — xPq{H)] exp (— isu;Fsu,), where x is the occupa- 
tion probability of |e > at the origin of the switching 
plateau 14J. Depending on the process, the state of the 
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FIG. 4: (a) Sequence of two successive readouts used 
to analyze the qubit state after the first measurement, 
(b) Conditional probabiUty for the measurement outcome 
L after a switching event in the first readout. (c) 
P^{H),Po{H),Po{e,H) and (d) oscillator switching rates 
TsiD , rf tu and qubit relaxation rates Vf , Vf as a function of 
the switching plateau driving amplitude Isw 

qubit after switching is either \g > or |e >. Therefore, 
it is possible to discriminate between the two processes 
by measuring the qubit state after a switching event. At 
the end of the switching plateau P(e, H) is given by: 

+ [l_^pe(^)]_^[,-*.„r.„ _,-*.„rf] (2) 

If relaxation is the only process which can lead to os- 
cillator switching then P{e, H) — 0. In equation ([2]), the 
first term proportional to PQ{e, H) corresponds to switch- 
ing events during regime I, whereas the second term pro- 
portional to corresponds to switching events during 
regime II. 

The occupation probability of |e > after a switching 
event, given by the conditional probability P{e\H) — 
P{e, H) / P{H), is measured with a second, successive 
readout. The conditional probability P{L\H) for a mea- 
surement outcome L after a switching event is P{L\H) = 
AP{e\H) + B, where A = 0.9 exp {-t hoi d^f) includes 
the contrast of the second readout and the qubit relax- 
ation, and B — 3.5 10~^ is the finite measurement error 
when the qubit is in \g >. When the oscillator is in 
the H state, qubit relaxation during the holding plateau 
is very strong. Therefore, the remaining signal is very 
small P{L\H) w lO^'^. A key point for our analysis is 
the very high contrast of the readout allowing a resolu- 
tion of P{L\H) as small as 10~^. T^^ and PQ{e,H) are 



extracted from a fit of P(e, H) to equation ([2|. The decay 
rates Tgw and are fixed (extracted previously), and 
P§{e,H) and T^^ are fitting parameters. The obtained 
values of Pq (e, H), T^^ and are shown in figures |4j:,d. 
For comparison, Pq{H), P§{H), Ff^, and F^^ obtained 
previously are shown as well. 

The ratio F^^^/Ff^, determines the intrinsic readout 
fidelity. The solid black lines in figure |4|i are fits of 
Fj^ and Ff.^, using the oscillator escape rate equation 
as given in [TT], from which we extract Ib2{g) = 108 nA 
and /f,2(e) = 136 nA. Extrapolating F^^, the intrinsic 
readout fidelity would be 99.9 %. However, in practice 
we observe a saturation of F^^^ at approximately 4 MHz. 

The qubit relaxation rate F^' can directly limit the 
readout fidelity. The independent determination of 
and F^^ is an important result of this paper. At the 
driving amplitude /°^* used for the measurement, F^' is 
higher than F^^. Therefore, for measurement times be- 
yond the oscillator transient period, qubit relaxation is 
the main process limiting the readout fidelity. 

To summarize, we have introduced a simple model 
to characterize the qubit relaxation and the oscillator 
switching process during measurement of a flux qubit. 
Qubit relaxation increases significantly with increasing 
driving strength with the oscillator in the low-amplitude 
state; it jumps to a much higher rate when the oscilla- 
tor switches to its high- amplitude state. Corrected for 
qubit relaxation and initial qubit preparation errors, the 
readout fidelity is at least 95 %. 

We thank F. Wilhelm, I. §erban, S. Ashhab, Y. 
Nazarov and H. Wei for useful discussions. This work 
was supported by the Dutch Organization for Fundamen- 
tal Research on Matter (FOM), E.U. EuroSQIP, the E.U. 
Marie Curie program, and the NanoNed program. 



[1] M. Devoret, A. Wallraff, and J. Martinis, cond- 
mat/0411174. 

[2] A. Lupa§cu, et ai, Phys. Rev. Lett. 96, 127003 (2006) 
[3] A. Lupa§cu, et al, Nature physics 3, 119 (2007) 
[4] N. Boulant, et al, Phys. Rev. B 76, 014525 (2007) 
[5] N. Katz, et al, Science 312, 1498 (2006) 
[6] A.N. Korotkov, Phys. Rev. B 60, 5737 (1999) 
[7] E. Il'ichev, et ai, Phys. Rev. Lett. 91, 097906 (2003) 
[8] A. Lupa§cu, et al, Phys. Rev. Lett. 93, 177006 (2004) 
[9] J.E. Mooij, et al, Science 285, 1036 (1999) 
[10] L Siddiqi, et al, Phys. Rev. Lett. 93, 207002 (2004) 
[11] M.L Dykman, Physica 104A, 480 (1980) 
[12] L Serban, PhD thesis Ludwig-Maximilians University, 

Munich, Germany, 2008 
[13] Taking into account the period T — 3.2 ns and the decay 
time r = 15 ns of the Rabi oscillations, we estimate the 
preparation fidelity of |e > to be 95%. 
[14] Due to the initial preparation errors 5% and the qubit 
relaxation during the rising part of the readout pulse 5%, 
X = 0.90. 



